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We consider the Grand Unified SU(5) model with a small or vanishing cubic term in the adjoint 
scalar field in the potential. This gives the model an approximate or exact Z2 symmetry whose 
breaking leads to domain walls. The simplest domain wall has the structure of a kink across which 
the Higgs field changes sign ($ — $) and inside which the full SU(5) is restored. The kink is 
shown to be perturbatively unstable for all parameters. We then construct a domain wall solution 
that is lighter than the kink and show it to be perturbatively stable for a range of parameters. 
The symmetry in the core of this domain wall is smaller than that outside. The interactions of 
the domain wall with magnetic monopole is discussed and it is shown that magnetic monopoles 
with certain internal space orientations relative to the wall pass through the domain wall. Magnetic 
monopoles in other relative internal space orientations are likely to be swept away on collision with 
the domain walls, suggesting a scenario where the domain walls might act like optical polarization 
filters, allowing certain monopole "polarizations" to pass through but not others. As SU(5) domain 
walls will also be formed at small values of the cubic coupling, this leads to a very complicated 
picture of the evolution of defects after the Grand Unified phase transition. 



I. INTRODUCTION 

Topological defects can be produced at a symmetry 
breaking phase transition and would be long-lived relics 
of the symmetric phase. If topological defects were pro- 
duced during a phase transition in the very early uni- 
verse, they could survive until the present epoch and thus 
provide a window to the very early universe. The lack 
of observable defects in the present universe helps place 
strong constraints on particle physics model building and 
early universe cosmology. 

A prototype symmetry breaking relevant for Grand 
Unified particle physics is 

SU(5) -> [SU(3) X SU(2) X \J{1)]/Ze . 

The corresponding phase transition would produce mag- 
netic monopoles. If the only factors affecting the evo- 
lution of the monopoles are the sub-luminal expansion 
of the universe and monopole-antimonopole Coulombic 
interactions, the monopole abundance grossly violates 
the observed absence of monopoles in the present uni- 
verse. The monopole over-abundance problem is solved 
by invoking superluminal universal expansion [i.e. infla- 
tion ^) or by extending the particle physics model so 
that the U(l) symmetry gets broken for a short dura- 
tion leading to confining forces between monopoles and 
antimonopoles |0 and thus enhancing their annihilation 
rat(^ Recently we have investigated the possibility 



* There is another possibility - that the Grand Unified 
phase transition never occurred and hence there never was 
a monopole over-abundance problem m . 



that the Grand Unified phase transition may also have 
produced a network of domain walls together with the 
magnetic monopoles. These walls would interact with the 
monopoles and sweep them away, reducing their abun- 
dance to an acceptably low level. It is to pursue this 
scenario in greater detail that we now study the struc- 
ture of domain walls in the SU(5) model. 
The bosonic sector of the SU(5) model is: 



where $ is an adjoint of SU(5), D^^^ 

are the gauge fields, and the potential is 



(1) 

;iven by: 



= -m2Tr($2 



/i[Tr($2)]2 

ATr($4)+^Tr($3)-l^o , (2) 



where, Vq is a constant that we will choose below. The 
SU(5) symmetry is broken to [SU{i) x SU[2) x U{l)]/Z(i 
if the Higgs acquires a vacuum expectation value (VEV) 
equal to 
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diag(2,2,2,-3,-3) 



where 



rj = 



m 

7^ 



\' = h 



(3) 



(4) 



(5) 



For the potential to have its global minimum at $ = (f>o, 
the parameters are constrained to satisfy: 

A > , A' > . (6) 
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For the global minimum to have V{^o) = 0, in eq. 
we set 

rr 4 



= tanh(o-z)$o 



(8) 



(7) 



The model in eq. (0) does not have any topological 
domain walls because the vacua related by $ ^ — $ are 
not degenerate. However if 7 is small, there are walls 
connecting the two kinds of vacua that are almost topo- 
logical. In our analysis we will set 7 = 0, in which case 
the symmetry of the model is SU(5)xZ2 and an expecta- 
tion of $ breaks the Z2 symmetry leading to topological 
domain walls in addition to the magnetic monopoles aris- 
ing from the SU(5) breaking. 

In this paper we will study the domain walls present 
in the SU(5)xZ2 model. The simplest kind of domain 
wall is the kink that has been studied in a single scalar 
field model with Z2 -> 1 (e^. [§). In |] we studied the 
interaction of the SU(5) kink with magnetic monopoles 
and found that the monopoles spread out along the kink 
on collision and never pass through. This confirmed the 
conjecture in Ref. [Q that kinks could sweep away mag- 
netic monopoles. However, the investigations of this pa- 
per show that the kink solution of the SU(5) x Z2 model is 
unstable to perturbations. The model contains another 
domain wall solution that is lighter than the kink and is 
perturbatively stable. The adjoint field does not vanish 
in the core of these new domain wall solutions and hence 
only a subgroup of the SU(5) is restored at the center. 
For this reason, the interactions of these domain walls 
with magnetic monopoles is expected to be much more 
complex (as compared to the kink), depending on the 
particular group orientation of the monopole relative to 
the wall. 

We will begin our analysis in Sec. || by constructing 
the kink and performing the stability analysis. Then in 
Sec. [II we will proceed to construct the domain wall 



in the model, prove that it is lighter than the kink, and 
show that it is perturbatively stable for a range of pa- 
rameters. In Sec. IV we will consider the interaction of 



monopoles and domain walls and show that a monopole 
whose orientation in group space is aligned with a col- 
liding domain wall, will pass through and not get swept 
away. We further conjecture that monopoles that are 
misaligned with the domain wall will be swept away but 
have not yet shown this. We draw an analogy of the 
sweeping out process with that of a polarization filter 
that "sweeps out" orthogonally polarized light and only 
lets through a certain polarization. 



II. SU(5) KINK: SOLUTION AND STABILITY 

The kink solution is the Z2 kink along the $0 direction 
(see eq. (||)). Therefore: 



with a = ni/y/2 (see eq. (||)), and all the gauge fields 
vanish. It is straightforward to check that solves the 
equations of motion with the boundary conditions ^(z = 
±00) = ±$0- 

As is well-known the mass (per unit area) of the 
kink is: 



Ah 



(9) 



Here we will examine the stability of the kink under 
general perturbations. So we write: 



(10) 



Since the kink solution is invariant under translations 
and rotations in the a;?/— plane, it is easy to show that the 
perturbations that might cause an instability arise from 
perturbations of the scalar field and can only depend on 
z. Therefore we may set the gauge fields to zero and take 
^ = '^{t,z). 

The Z2 kink is stable and hence we can restrict the 
scalar perturbations to be orthogonal to Further- 
more, since the stability of the kink to diagonal pertur- 
bations has already been studied in Ref. Q , we only have 
to consider perturbations that cannot be diagonalized by 
a global SU(3)xSU(2)xU(l) transformation that leaves 
the kink invariant. Therefore we can write: 
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(11) 



where are all generators of SU(5) that do not commute 
with $0. 

Next we analyze the linearized Schrodinger equation 
for small excitations tp^ = ^Q(z)exp{—iujt) in the back- 
ground of the kink: 



-d',~m^ + cl,Uz){h + XnWo^Lo^ro 



(12) 



where (jjk = tanh((Tz) and = 7/30. Since this equation 
is identical for excitations along any of the 12 directions, 
we can drop the index i. The kink is unstable if there is 
a solution to eq. (p^ ) with a negative w^. Substituting 
eq. (|) into eq. (1121) yields: 



+ m2(tanh^(crz) - l)]?/'o = ^^V'o 



(13) 



This equation has a bound state solution ipo oc sech.(az) 
with the eigenvalue = — to^/2. Since this result is 
independent of the parameters in the potential, we con- 
clude that the kink in SU(5) is always unstable. 



2 



III. DOMAIN WALL 

The domain wall solution is obtained if we choose the 
gauge fields to vanish at infinity and the scalar field to 
satisfy the boundary conditions: 



$(z = -oo) = 



V 



2V15 
5 V 



diag(2,-3,2,2,-3) 



V\lj^{h + T3) + j{Y-V5Xs) (14) 



and 



<^{z = +00) = $+ = ;^^diag(3, -2, -2, 3, -2) 



2V15 

= V^ih+r3)-^{Y-V5Xs) . (15) 

Here A3, As, T3 and Y are the diagonal generators of 
SU(5): 



A3 = idiag(l,-l,0,0,0) 



diag(l, 1,-2, 0,0) 



2^3 



(16) 



(17) 



+ h{a'' + b^)]c (23) 



d" = [-m^ + {h + ^)d' + {h+^)ia 
30 5 



2 , 52^ 



9A, A6 , 2 &^ 



(24) 



where primes refer to derivatives with respect to z. For 
reference, the kink solution (eq. (||)) corresponds to 
a{z) = = b{z) — c{z) and d{z) = ?7tanh(crz). 

The equations of motion for b and c and can be solved 
quite easily: 



b{z) = —V5d{z) , c(z) = a{z) 



(25) 



This is consistent with the boundary conditions in eqs. 
(|l|) and dH). In addition, we require 

a{z = ±00) = +^7^^ , d{z = ±cx)) = t| • (26) 
Then the remaining equations we need to solve are: 



6h+^X]d^ 



a (27) 



T3 = -diag(0, 0,0, 1,-1) , 



(18) 



d" 



-m'+[ 6h+^X^d^+(2h+^]a^ 



d (28) 



Y 



1 
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diag(2,2,2,-3,-3) . 



(19) 



These equations can be written in a cleaner form by 
rescaling: 



Note that a local SU(5) transformations can be used to 
rotate <!>+ into — $~ so that the boundary conditions are 
like those of the kink with <i>(z = +00) = —^{z = —00). 
However, then the solution for the domain wall will not 
be diagonal at all z. We prefer to use the above boundary 
conditions so that the solution is diagonal throughout. 

The domain wall solution can be written as 



A{z) = 
Then 

A" 



-- , L»(z) = -6- , Z = mz. (29) 
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-1 + ^l^D^ + (i±^A^' 
5 5 



A 



(30) 



^Dwiz) = a{z)X3 + biz)Xs + c(z)t3 + d{z)Y . (20) 

The functions a, b, c, and d must satisfy the static equa- 
tions of motion: 

a" = l-m^ + {h+ ^)d^ + {h + ^){a^ + b") + h^]a 
z 



b" ^ 



2X 



2Xabd 

~7r 



A. 



(21) 



+ {h+ —)d' + {h + -){a' + b') + hc']b 



Xd , , 



(22) 



D" 



-l+pL'2 + (l-p)yl2 



(31) 



where primes on A and D denote differentiation with 
respect to Z, and 



V 



1 



5A 



(32) 



Note that p G [1/6, 00) because of the constraints in eq. 
(^). The boundary conditions now are: 



A(z = ±oo) = +l, D(z = ±oo) = ±l 



(33) 



This system of equations has been solved by numer- 
ical relaxation and a sample solution is shown in Fig. 
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To find an approximate analytical solution, assume 
that |j4"/A| << 1 is small everywhere. This assumption 
will be true for a certain range of the parameter p which 
we can later determine. Then the square bracket on the 
right-hand side of eq. (pOj) is very small. This gives: 



A: 



1 - 



1/2 



We insert this expression for A in eq. 
the kink-type differential equation: 



(34) 



and obtain 



where 



D" = q[-l 



&P-1 



6A 
A + 60A' 



and the solution is: 



D{Z) ~ tanh( 



(35) 



(36) 



(37) 



The parameter q lies in the interval [0, 6]. For q = I {i.e. 
p = 1) it is easy to check that this analytical solution is 
exact. 

We can now check that our assumption << 1 is 

self-consistent provided p is not much larger than a few. 

The energy density for the fields A and D can be found 
from the Lagrangian in eq. (|l]) together with the ansatz 
in eq. (pO[), the solution for b and c in eq. ( |25| ) and the 
rescalings in eq. (p9[). The resulting expression for the 
energy per unit area of the domain wall is: 



M 



DW 



12A' 



dZ[5A'^ + D'^ + ViA,D)] (38) 



where, 



ViA, D) = -5^2 -D^+ (^±11^4 



^D^ + {l-p)A^D'' +?, . (39) 



The energy can be found numerically. However, here we 
will find an approximate analytic result. We can insert 
the approximate solution given above in eq. ( ^ ) but this 
leads to an expression that is not transparent. Instead it 
is more useful to consider another approximation for A 
and D : 



A ~ 1 , £> ~ tanh^^|z 



(40) 



(This approximation is exact for p = 1.) A straightfor- 
ward evaluation then gives: 



MoWapprox = ^^^^^ 



(41) 



where, is given in eq. (^. 

We can now compare the domain wall energy to the 
kink energy. If the domain wall is the least energy so- 
lution for the given boundary conditions, the energy of 
the exact solution for the domain wall will be bounded 
above by the energy of the approximate solution. Note 
that this will be true even if the approximation used to 
find the analytical solution is not good. Hence this sim- 
ple argument shows that the domain wall is lighter than 
the kink for p < 36, or for h/X > —6.94/30. A numerical 
analysis shows that the domain wall is lighter than the 
kink even in the range -6.94/30 > h/X > -7/30. There- 
fore the domain wall is lighter than the kink for the full 
range of parameters specified in eq. (||). 



FIG. 1. The domain wall solution for A = 1 and h = —0.2 
(p — 2.25). The solid line shows a{z) and the dashed hne 
shows d{z). 

Next we study the stability of the domain wall solu- 
tion. It is easy to show that the solution is stable to 
diagonal perturbations, so here we focus on off-diagonal 
perturbations. We write: 



(42) 



where N"' are the non-diagonal generators of SU (5) and 
■0" are small perturbations satisfying the boundary con- 
ditions tp°'{±oo) = 0. Let us first consider the contribu- 
tion to the energy density due to fields tp°'. To second 
order in perturbations the contributions from different 
modes, labeled by index a, do not couple. A more de- 
tailed analysis shows that the mode corresponding to 



iVi = 



^0 1 0\ 

1 0' 





Vo o/ 



(43) 
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is one of the 8 modes that are most unstable. Let ij) 
be any of these 8 fields. The contribution to the energy 
density due to ijj is: 



1 



A 
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{a + —d) + higher order terms 



(44) 



where a and d are defined by eq. (|2(]|). As in the case of 
the kink, we are interested in the linearized Schrodinger 
equation for small excitations "0 = 'ipo{z)^xp{—iujt) in the 
background of the diagonal domain wall solution. Eq. 
(H) leads to 



(Gh 



l.X)d' + (2h+^)a^]i,o 



Comparing this with eq. (|2^) allows us to write: 
a" 

- "00 H ■00 = w^Vo ■ 



(45) 



(46) 



If a" /a = 0, as happens when p ~ 1, then there are no 
non-trivial solutions to eq. (^) that satisfy the correct 
boundary conditions. Therefore, the diagonal domain 
wall solution is stable for at least one choice of parameters 
in the potential, namely, for p = I. By continuity there is 
a range of parameters around p — 1 for which the domain 
wall is perturbatively stable. 



IV. INTERACTION WITH MONOPOLES AND 
DISCUSSION 

To understand the interaction of the domain wall with 
magnetic monopoles, it is first useful to understand the 
structure of the domain wall core. Since a{z) is non-zero 
inside the domain wall, ^Dwiz = 0) = a(0)(A3 -I- ra) cx 
diag(l, — 1, 0, 1,-1). Therefore, the symmetry inside the 
core is K = SU(2)xSU(2)xU(l)xU(l). The first 
SU(2) factor arises due to rotations in the 2x2 block 
with the entries equal to 1 in ^dw (first and fourth rows 
and columns). The second SU(2) factor is due to the 
block with entries equal to -1 (second and fifth rows and 
columns). The two U(l) factors arise since there are two 
diagonal generators of SU(5) aside from those already 
accounted for in the two SU(2) factors, that commute 
with ^Dwiz — 0). (We are ignoring any discrete factors 
that might be present.) The symmetry group K within 
the domain wall is to be contrasted with the full SU(5) 
symmetry which is restored within the kink. The fact 
that only a subgroup of the SU(5) symmetry is restored 
in the core of the wall means that the interaction of the 
monopole will now depend on the particular embedding 
of the monopole in SU(5) and its orientation in internal 
space relative to the domain wall. 



Consider a magnetic monopole whose winding lies in 
the fourth and fifth rows and columns of $. Staying close 
to the notation of Q we write the scalar field of such a 
monopole as: 



*M(r) = P(r) J2 ^"^a + M{r)\'^ + N{r)Y , (47) 

a=l 

where {to} are the SU(2) generators (see eq. (p^)) and 
A' 



^diag(I,-2,l,0,0) = ^A3 



2 ^ 



The non-zero gauge fields are: 
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X? 



a=l 



e%- — (l-if(r)) 



(48) 



The monopole profile functions, P{r), M{r), N{r) and 
K{r), are solutions of the static equations of motion with 
boundary conditions: 



/5 \/5 
— , M(oo) = rj— , 

iV(^) = I , K{^) = I. 
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(49) 
(50) 



When the monopole and the wall are very far from each 
other, the combined field configuration can be described 
by the following ansatze: 



M+DW 



+M{r) 



c(-cxo) ^ 

^ ' a— 1 

VS a{z') 
2 a(— CX)) 



N{r) 



d{z') 



d{—oo) 

1 b[z') 
2b{-oo) 



-Y 



(51) 



where z' — z zq and zq is the initial position of the 
domain wall. When r is small ^m+dw (eq. (|4^)) 

and when z' is small ^m+dw ^dw (eq- (20)) along 
the z-direction. The gauge fields are the same as for the 
monopole alone. We have purposely chosen the embed- 
ding of the monopole so that all interesting dynamics of 
the monopole-wall interaction is restricted to the fourth 
and fifth rows and columns of ^m+dw- This follows from 
the equations of motion and the commutation properties 
of the generators appearing in the ansatze (pT]). Let us 
then only consider the relevant part of the $ m+ dw ma- 
trix: 



2 a(-oo) \x + iy -z 



d{z') /I 
2VT5^''''rf(-oo) \Q 1 



zN{r 



(52) 
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The form of $2x2 suggests that the only field that is 
going to be considerably affected by the domain wall is 
N{r) because a is roughly constant in space. There is 
no angular dependence in the term with N{r) in eq. ( ^ ) 
and hence N{r) does not contribute to the winding of the 
monopole. Therefore, we do not expect the wall to affect 
the winding. Essentially the reason is that the SU(2) 
subgroup in which the monopole winding is located is not 
restored on the wall. We have checked that the monopole 
passes right through the wall explicitly in this case by 
numerically colliding the monopole and the wall. 

If the magnetic monopole winding lies in the first and 
fourth block, it will experience a region of restored SU(2) 
symmetry inside the domain wall and hence we conjec- 
ture that such monopoles will unwind on the wall. If 
our conjecture is correct, the domain walls behave sim- 
ilarly to optical polarization filters, allowing monopoles 
with certain internal space polarizations to pass through 
and annihilating other polarizations. The detailed anal- 
ysis of all possible monopole embeddings is a challenging 
project, both numerically and analytically, since one can- 
not avoid dealing with a large number of the fields present 
in SU(5). 

There is another possibility that is worth pointing out. 
If a domain wall and a magnetic monopole are misaligned 
in internal space, it may not be possible to superpose the 
two solutions so as to get a monopole and a domain wall 
together. (Such a situation is known to occur when at- 
tempting to construct multimonopole or multistring so- 
lutions.) Then it is likely that there will be a long range 
force between the domain wall and misaligned monopole 
that will bring them together. On coming together the 
monopole could get annihilated on the wall, or else in 
some cases, it may get aligned and then pass through the 
wall. 

Our considerations point to a very complicated af- 
termath of the GUT phase transition. Domain walls 
and magnetic monopoles would both be produced and 
would start interacting. The outcome of an interaction 
would depend on the internal space orientations of the 
monopole relative to the domain wall. Any given domain 
wall would be transparent to some monopoles but not to 
others. The relaxation of the system would depend on 
whether a monopole encounters a sufhcient number of 
randomly oriented (in internal space) domain walls, at 
least one of which might sweep it away. It remains to 
be seen if domain walls can provide a means to solve the 
cosmological monopole over-abundance problem. 
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